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ABSTRACT 

The final merger of a pair of massive black holes in a galactic nucleus is compelled by gravitational radiation. 
Gravitational waves from the mergers of black holes of masses (10 5 - 10 7 )(1 + z) _1 M at redshifts of 1-20 
will be readily detectable by the Laser Interferometer Space Antenna, but an electromagnetic afterglow would 
be helpful in pinpointing the source and its redshift. Long before the merger, the binary "hollows out" any 
surrounding gas and shrinks slowly compared to the viscous timescale of a circumbinary disk. The inner gas 
disk is truncated at the radius where gravitational torque from the binary balances the viscous torque, and 
accretion onto the black holes is diminished. Initially, the inner truncation radius is able to follow the shrinking 
binary inward. But eventually the gravitational radiation timescale becomes shorter than the viscous timescale 
in the disk, leading to a merged black hole surrounded by a hollow disk of gas. We show that the subsequent 
viscous evolution of the hollow, radiation pressure-dominated disk will create an ~ 10 43 5 (M/10 6 M Q )ergs s" 1 
X-ray source on a timescale ~ 7(1 +z)(M/10 6 M Q ) L32 yr. This justifies follow-up monitoring of gravitational 
wave events with next-generation X-ray observatories. Analysis of the detailed light curve of these afterglows 
will yield new insights into the subtle physics of accretion onto massive black holes. 
Subject headings: accretion, accretion disks — black hole physics — quasars: general — X-rays: galaxies 



1. INTRODUCTION 

Evidence is mounting that galactic spheroids contain mas- 
sive black holes (MBHs) in their nuclei. When two galaxies 
merge, their MBHs form a binary in the nu cleus of the new 
galaxy (Begelman, Blandford, & Rees 1980). The binary in- 
teracts with its stellar and gaseous environment and with other 
MBHs that collect in the same nucleus in multiple mergers. 
These interactions torque the binary and extract its orbital an- 
gular momentum and energy. The binary can thus be rendered 
so compact that gravitational radiation (GR) carries away the 
remaining orbital energy, inducing coalescence. For binaries 
with masses M = M\ +M2 < 10 7 M Q , the gravitational sling- 
shot ejection of stars is sufficient to guarantee coalescence in a 
Hubble time dMilosavlievic & Merrittll2"b 03). GR emitted by 
these binaries shortly before and during coalescence will be 
detected by the space-based gravitational wave detector Laser 
Interferometer Space Antenna (LISA). 1 The coalescence of 
isolated black holes in general relativity is not accompanied 
by observable electromagnetic (EM) emission. We here show 
that circumbinary gas can lead to a delayed EM afterglow. 

Of interest to LISA are mergers of binaries with masses 
M < 10 7 M Q of arbitrary mass ratio q < 1. There is abun- 
dant evidence for dense gas in galactic nuclei. Geomet- 
rically thin molecular disks have bee n detected in water 
maser emission in Se yfert galaxies (e.g., Mivoshi et al.ll995l 
iGreenhill et al.ll2003l) . The Galactic nucleus contains a 4 x 
10 6 M(7) M BH surrounded by an ~ 10 4 M Q molecular gas torus 
(e.g., Qackson et al. 1993). Massive accretion disks must be 
present in quasars and narrow-line Seyfert I nuclei to explain 
their luminosities. 

In general, the specific angular momentum of inflow- 
ing gas exceeds that of the binary, while the gas temper- 
ature is below the virial temperature GM/j,m p /2ak, where 
\i is the mean molecular weight in units of the proton 
mass OT p , a is the binary's semimajor axis, and k is the 
Boltzmann constant. The sub-virial gas settles into a ro- 

1 http://lisa.jpl.nasa.gov. 



tationally supported, geometrically thin circumbinary disk. 
If the disk is inclined with respect to the binary's orbital 
plane, the quadrupole component of the binary's gravita- 
tional potential causes differ ential precession and the warp- 
ing of the disk. As in the Bardeen & Petterson ( 1975) mech- 
anism, the warp dissipates viscously, resulting in a disk 
in the binary's orbital plane tLarwood & P apaloizoul 119971 
Ivanov. Papaloizou. & Polnare vl1999h . 

The disk is truncated at an inner edge where 
gravitatio nal torques and vis cous stresses balance 
(e.g., lArtvmowicz & Lubowl 119941) . Surface density in 
the h ollowed region is much smalle r than in the disk 
(e.g., iGunther. Schafer. & Klevl 12004 . While material 
does peel off the inner edg e and flow across the hole 
( Artvmowicz & Lubow 1996), the accretion rate for mod- 
erate mass ratios (q > 0.01) is only a fraction (< 10%) of 
the rate at which the disk would be accreting without being 
torquedby the binary (Lubow. Seibert. & Artvmowicz 1999). 
As the binary's semimajor axis decays owing to stellar or 
gas dynamical processes (including angular momentum 
extraction by the circumbinary disk itself), the disk's inner 
edge spreads inward viscously while maintaining an approxi- 
mately constant ratio of inner edge radius to semimajor axis, 
fedge/a ~ 2. In the final stages of in-spiral, however, the time 
for decay of the semimajor axis is set by GR and decreases 
rapidly with a (equation^. When this becomes shorter than 
the (viscous) time for the inner edge of the disk to spread 
inward, the binary conti nues toward merger while th e disk 
structure remains frozen ( Armitage & Nataraian 2002). This 
was also noted by iLiu. Wu. & Cad d2003h . who suggested 
that the formation of an inner hole in the disk preceding 
MBH coalescence causes the interruption of jet activity in 
double-double radio galaxies. 

Vertical support in the central parts of MBH accretion disks 
is dominated by radiation pressur e, while opacity is domi- 
nated by electron scattering. The Shak ura & Sunvaevl (119731) 
q-disk s in this regime are thermally and viscously unstable 
(Lightman & Eardlev3 Q974l) if viscous stresses scale with to- 
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tal pressure. However, they are stable if the stre sses sc ale with 
the gas pressure alo ne ISa kimofo & Coronitilll981l) . Such 
a scaling can occur (Turner 2004) if the distance photons 
diffuse per orbit is about equal to the scale of magnetic 
fields produced b y the magnetorotational instability (MRI; 
Balbus & Hawle*vl ri991l) that facilitates angular momentum 
transport. The radiation then partly decouples fr om the tur- 
bulence driven by the MRI (iTurner et al. 1 120031) . These 
disks are cl umpy and porous to the radiation produced within 
SB e gelmanl2002l and references therein). 

We here describe the observable signatures of such cir- 
cumbinary disks in MBHs. In §|5J we present a scenario for 
the combined evolution of an MBH binary and its circumbi- 
nary accretion disk immediately preceding coalescence. In 
§|3] we discuss observable signatures. 

2. A MODEL FOR BINARY-DISK EVOLUTION 

In the early stages of the GR-driven orbital evolution, the 
binary orbit is circularized by GR. Furthermore, unless q <C 1, 
the timescale ( Pete rs! 1 9641) 

t a - 5 cV (l+q)2 m 

gr_ da/dt 64 G 3 M 3 q 

on which the binary's semimajor axis a decays because the 
emission of GR greatly exceeds the viscous timescale ? v j sc = 
(2/3)r 2 /^(r) for r ~ a. For realistic disks (see below), the 
effective kinematic viscosity v is a weak function of radius, so 
? v isc oc r 2 approximately. Thus, any gas inside the binary orbit 
should long ago have been accreted (for r <C a) or expelled 
(for r ~ a), leaving the MBHs in a gas-free "donut hole." 2 

Gas outside the binary will attempt to accrete but be pre- 
vented by torques from the binary (§Q. The inner edge of the 
disk lies at r e d ge = 2Xa, where A > 0.5 is a parameter of order 
unity. In this regime, the vis cous torque in the disk 3-7rr 2 z/£f2 
is constant with radius (e.g.,|Pringle 1991). Here £ is the sur- 
face density and Q = (GM /r 3 ) 1 ! 2 is the local angular velocity. 
As the binary shrinks, the inner parts of the circumbinary disk 
pass through a time sequence of nonaccreting, constant torque 
configurations. In the outer parts of the disk it can occur that 
?gi < Wisdr), where r is a radius in the disk. Then material 
outside r loses viscous causal connection with the inner edge 
and decouples from the binary torque. 

Since f gr oc a 4 , while the viscous timescale at r e d ge approx- 
imately scales as f v i sc oc a 2 , the shrinking binary/disk system 
eventually reaches a state in which f gl < f v isc('"ed g e) I ie., the 
whole disk has decoupled from the rapidly shrinking binary. 
The inner parts of the disk begin to evolve as a standard, 
zero-torque disk. When the inner edge decouples, the ra- 
dial profile closely resembles the non-accreting disk with a 
sharp edge. The inner edge reaches the center in time f s h ~ 
/3?visc(red ge ), where (3 = max{0.1, [dlnS/c/lnKredge)] -1 }. The 
factor 0.1 is appropriate for an infinitely sharp disk 
dLvnden-Bell & Pringldll974T) . To estimate f s h and the self- 
consistent a for which f gl ~ f s h, we must determine the in- 
ternal structure and viscosity of the inner disk at the time of 
decoupling, to which problem we turn now. 

In the stable a-model, the kinematic viscosity is given by 
v ~ |a ga s-P g as/ where f gas = pkT /pm p is the gas pressure, 
and p is the density. To prevent confusion, we define a gas , 
related to the usual definition via a g a S ^gas = ccftotai, where 

2 This is in contrast to the q <§C 1 case considered by Armitaae & Nataraian 
1 2002), who assumed that there was still gas accreting onto the black holes. 



^totai = Pm&+Pgas + Pmag is the sum of the radiation, gas, and 
magne tic pressures . In a simulation of one such accretion 
disk, ITurner! ( 120041) measures a = 0.0013 and P nd /P g ^ = 14 
while P mag <C P m d, implying that a gas ~ 0.02 in his simula- 
tion. We defer modeling the detailed vertical structure of the 
disk to a follow-up paper; relations presented below are ap- 
proximate. 

The rate of local dissipation per unit area of the disk 
equals Q ~ (9/4)^EJ7 2 . When we ignore horizontal advec- 
tion, the flux F = Q/2 must be emitted from each surface. 
If the surface is emitting as a blackbody, the flux equals 
Fbb = 4oT 4 /3t, where T is the midplane temperature, r is 
the total optical depth, and a is the Stefan-Boltzmann con- 
stant. The spectrum differs from that of a blackbody since 
the photons at different frequencies are thermalized at differ- 
ent depths. The modified, "graybody" spectrum has emergent 

flux F v ~ 7rei/ 2 (l + e l J 2 )~ y B v , where B w is the Planck func- 
tion and t v = K a bs,i//(K a bs,i/ + Kes) is the ratio of the absorption 
to the total opacity, and K es is the electron scattering opacity. 
The quantities B v and e„ are evaluated at the bottom of the 
thermalization photosphere (TP); we denote the temperature 
and the density there by T v and p„, respectively. 

The scale height of the photosphere is given by h v ~ 
["fPmd(Tv)/ p v ] 1 / 2 /Q, where 7 w 4/3 is the adiabatic index. 
We set the absorption opacity equal to the (mainly bound- 
free) opacity k„r;3x lO^gwpT' 1 l 2 £,~ 3 {\-e~^) cm 2 g -1 , with 
p and T in cgs, where £ = hv/kT, h is the Planck constant, and 
ghi ~ 1 is the gaunt factor (k„ is the free-free opacity scaled 
up by the ratio of the bound-free to the free-free Rosseland 
mean opacities at solar metallicity). By definition, the effec- 
tive optical depth ~ (Kabs,!/« es ) Pv n v = 1 at the bottom of TP. 
This is solved for p„, which we substitute back in K a bs to find 

that e 2 = 1.4 x 1 42 ^ b f ^ 2 7l7 15/2 ^" 3 ( 1 -e*). To estimate the 
degree to which blanketing by the TP modifies the integrated 
emitted flux, we evaluate e„ at the Wien frequency (£ « 2.8). 
We relate T v to the midplane temperature via T v = t -1 ' 4 7\ 
where r is the electron scattering optical depth between the 
photosphere and the midplane. This yields 

F gb = / F v dv — — ~ -vYSl , (2) 

Jo iT 1 + V e ° 

where e ~ 2.5 x 10 20 f2T 15//16 r" 15//4 and the last approximate 
equality in equation (0 follows from identifying f g b with half 
of the power generated in the disk. The disk edge has e < 1 at 
decoupling, resulting in a higher midplane temperature than 
in the blackbody disk. 

The optical depth is given by r = QnY,, where k is the 
opacity (electron scattering K es 1=3 0.4 cm 2 g" 1 , or Kramer's 
Kabs ~ 1-6 x 10 24 pr -7 / 2 g" 1 ) and 9 < 1 is a "porosity" cor- 
rection factor. We introduce 9 to account for the shortened 
effective optical depth in an inhomogeneous disk in which ra- 
diation escapes through low-density domains. The correction 
can also be used to approximate disks in which a fraction of 
the turbulent ma gnetic energy is dissipated in surface layers 
iMiller & $tone|2 p00l) o r in the presence of photon-bub ble in- 
stability (Gammie 1998). We estimate that 9 sa 0.2 in Turner 
(2004). 

The surface density is thus far unspecified. One (admit- 
tedly artificial) way to parametrize £ is to fix the accretion rate 
M = 3nvYi that the disk would have at its inner edge if it con- 
tained a single black hole instead of a binary (the disk is not in 
viscous equilibrium, and the accretion rate varies with radius). 
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The accretion rate can be expressed in units of the Eddington 
rate M e dd = 4-nGMm v /r)caj, where r\ ~ 0. 1 is the radiative ef- 
ficiency and <tt is the Thomson cross section; M e( jd should not 
be confused with the local Eddington limit at the disk edge. 
We define the dimensionless parameter m = M/M^. There 
is no a priori requirement for m < 1; in fact, following de- 
coupling the accretion rate will likely be the largest anywhere 
among black holes of the same mass. The surface density is 
thus parametrized via 

4 GMm v . 
via = -m. (3) 

3 r)co~T 

Another way to characterize the surface density is via the total 
mass of the disk Mdisk (see Table out to the radius where 
the temperature falls below 10 4 K (see below). 

To determine the disk properties at the time of decoupling, 
we equate f s h with f gl , which with the expression for the 
viscosity determines the midplane temperature near the disk 
edge in terms of the semimajor axis. This is then used in 
equations (|2j and (0 to solve for the semimajor axis a, the 
edge surface density S, and the edge midplane temperature 
T at decoupling. Self-consistent solutions are in the radi- 
ation pressure-dominated and electron scattering-dominated 
regime. The results are summarized in Tabled We have de- 
fined a gas = 0.1cm, 13 = 0.1/3-1, r\ = 0.177-!, M = 10 6 M 6 M Q , 
and 9 = 0.26>o.2. 

We calculate the disk scale height-to-radius ratio h/r = 
£ jrp. For the fiducial choice of parameters a~\ = (3-\ = rj-i = 
\ = m= M§ = q = #0.2 = 1 , the disk is marginally geometrically 
thick at the inner edge. We ignore horizontal and vertical ad- 
vection but do provide an estimate of the advected-to-radiated 
heat flux ratio Q a dv/Qrad- Ignoring radiative loss, we find 
the advected flux equals Q a d v = (3/4){4^ p - H^jfvYiP^/ pr 2 , 
where j p> t = dln(p, T)/d\nr\ in our disk, {4^f p - 127t) ~ 8. 
In the fiducial disk, horizontal advection is competitive with 
radiative diffusion, implying that the true disk is geometri- 
cally thinner than ours, as in the "slim-disk" solutions of 
lAbramowicz et all (Q988). Just how large could the param- 
eter m be? The disk must be thin at the edge (/i/r) e d ge -C 1, 
implying (we ignore the weak dependence on a, M, and p) 

m « 1.4/30fr ? - 1 A 1 - 56 [4^/(l + ^) 2 ]°- 39 ^ 49 ) (4) 

which allows for m 3> 1 when q <C 1. A trustworthy 
disk model should also be gravitationally stable, with Q = 
fi 2 l%Gp > 1 throughout. We have checked that the disk struc- 
tures for M < 10 6 M Q have Q > 1 inside the radius where T 
falls below 10 4 K and the disks become susceptible to ioniza- 
tion instability. We do not extrapolate our disks beyond this 
radius. 

3. DETECTION OF THE ELECTROMAGNETIC AFTERGLOW 

The pre-coalescent circumbinary disk is expected to be lu- 
minous in IR with a negligible X-ray counterpart (Figure 
Unfortunately, this IR source may be co nfused with the 
light of the host galaxy. Barth et al. (2004) recently stud- 
ied the dwarf galaxy POX 52, which appears to contain a 
1.6 x 10 5 M(7) black h ole accreting at about the Eddington rate 
iGreene & Hol20 04 list other similar candidates). The optical 
luminosity of the POX 52 black hole is about the same as that 
of the galaxy, so it seems unlikely that the MBH binary can 
be identified electromagnetic ally before the merger. 

The black holes merge f gr /4 after decoupling, where t $I is 
evaluated at decoupling. The gravitational wave losses during 



Table 1. Disk at Decoupling" 



Variable 


Factor 


Ot- 1 


f7_] /m 


X 


M 6 




0Q2 


a/(GM/c 2 ) 


117 


-0.34 


0.24 


0.70 


0.08 


0.42 


-0.08 


T (10 6 K) 


1.7 


0.19 


-0.86 


-1.95 


-0.28 


-0.49 


0.30 


fsh (yr) 


9.4 


-1.36 


0.98 


2.80 


1.32 


1.7,0.7 


-0.34 


h/r 


0.46 


0.76 


-2.43 


-3.80 


-0.12 


-0.95 


1.19 


Aad/^gas 


2600 


1.67 


-4.25 


-7.35 


-0.04 


-1.84 


2.17 


Kes /^abs 


46,000 


1.76 


-4.68 


-8.32 


-0.18 


-2.08 


2.32 


2(10 5 ) 


4.5 


2.12 


-2.41 


-6.61 


-1.44 


-1.65 


1.53 


Qadv/ Srad 


0.44 


1.52 


-4.86 


-7.60 


-0.24 


-1.90 


2.38 


e (lO" 3 ) 


1.4 


-0.84 


2.37 


4.20 


0.08 


1.05 


-0.21 


Mdisk («©) 


96 


-1.17 


-0.88 


-1.15 


2.04 


0.21 


-0.04 



a Variable in column 1, defined in the text, equals the factor in column 2 
times the product of column head parameters raised to the exponents indi- 
cated in columns 3-8, where we assume /i = 0.6. All quantities except for 
Mdjsk are evaluated at the edge of the disk. 
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FIG. 1 . — Representative thermal disk spectra before the coalescence (thin 
lines) and after (thick lines). The spectra are modified blackbody spectra of 
Eddington-limited (m =1) accretion disks around black holes of mass M = 
(10 4 , 10 5 , lO 6 )A/0 (from right to left) and dimensionless spin parameter a* = 
0.9 (we ignored Doppler broadening). Only the emission from rings in the 
disk with r < 1000GM/c 2 was take n into account . The se crude spectra are 
compatible with detailed models of Hubenv et al. 1 2001 ). We do not show 
thermal emission from any accretion disks around individual black holes that 
may be fed by the gas flow crossing the hollowed region. 



merger immediately perturb the potential in which the sur- 
rounding gas orbits, giving a (weak) prompt EM signature. 
Subsequently, the inner edge of the accretion disk migrates 
inward on timescale ~ f s h ~ ? gl , arriving at the merged black 
hole ~ |f s h after the merger. Its arrival is accompanied by 
rapid accretion and the activation of an X-ray active galactic 
nucleus (AGN). We consider these two distinct types of EM 
signatures in turn. 

The few percent reduction in the total black hole mass 
due to gravitational wave losses excites a weak axisymmet- 
ric wave in the disk. More interestingly, t he coalescence may 
be accompanied by radiation recoil (e.g., Merritt et al.l l2004t 
and references therein). The velocity of the recoiled black 
hole, v lec < 300 km s _I , is much smaller than the orbital ve- 
locity of the disk at the inner edge, ~2x 10 4 km s -1 . Disk 
material at radius r remains entirely bound to the black hole if 
its orbital velocity prior to the recoil (GM/r) 1//2 is larger than 
~ 2.4v rec , which is true for the inner disk. The recoil drives 
waves and warps in the outer disk. If these result in shocks or 
obscuration, they could result in an observable EM signature. 

More certain is the afterglow due to the viscous migration 
of the inner edge of the disk and consequent turn-on of X- 
ray emission from a rapidly accreting AGN. LISA will de- 
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M (M ) 

FIG. 2. — Redshifts at which emission from an Eddington-accreting (m = 1) 
black hole of mass M can be detected with 10 5 s exposure using the instru- 
ment ACIS-I on Chandra (dark gray) and the Hafnium superconducting tun- 
neling junction (STJ) detector on XEUS before mission upgrade (light-grey). 
We assume that 10 counts constitute detection and ignore absorption and con- 
fusion. The disk spectra were thermal modified blackbody spectra for a nearly 
maximally spinning black hole, a* = 0.9. We assume the standard ACDM 
cosmology with Hq = 70 km s _1 Mpc~', f2 m = 0.3, and £lj± = 0.7. 



tect with a high signal-to-noise ratio the coalescence of black 
holes with masses (10 5 - 10 7 )(1 +z)~ l M Q at redshifts z < 20. 
Can the EM afterglow also be detected? From Table \l\ we 
estimate the observed interval Af ~ |(1 + z)t s u between the 
two signals. For example, for the merger of two 10 5 M Q black 
holes at z = 5, assuming a.-\ = /3_i = ?y_i = m = Oq.2 = 1, we find 
Af ~ 2 yr for a graybody disk truncated at r e d ge = 2a. A strong 
recoil could shorten this. Most of the luminosity of the post- 
coalescent accretion disk around a (10 4 - 10 6 )M© black hole 
is emitted at rest-frame energies hv ~ (0.5-5) keV (see Fig- 
ureQ, which fall within the sensitivity windows of future X- 
ray detectors such as XEUS 3 and Generation-X. 4 XEUS could 

3 http://www.rssd. esa.int/index.php?project=XEUS. 

4 http://generation-x.gsfc.nasa.gov. 

5 http://www.stsci.edu/jwst. 



see this emission at z ~ 20 for > 10 5 M Q black holes. For 
smaller black holes, the maximum redshift for detection de- 
creases with the decreasing mass (see Figure|2j. 

At z < 20, the James Webb Space Telescope 5 is sensitive to 
rest-frame mid-UV to near-IR, which come from the disk radii 
exterior to the inner edge at decoupling (Figure[0. Therefore, 
the post-coalescence IR-optical flux will not be very different 
from the pre-coalescence flux. An exception would be if the 
merged black hole were enshrouded by a large column depth 
of gas and dust. Then, the UV-X-ray emission would be re- 
processed to longer wavelengths. Such a shrouded merger 
could be identified by a sudden increase in IR luminosity. 

We therefore propose that locations with detected GR 
signals from MBH coalescence be monitored in the 0.1- 
lOkeV band at high angular resolution. The angular res- 
olution of LISA ranges from several arcminutes to sev- 
eral degrees depe nding on the black h ole mass, mass ra- 
tio, and redshift (ICutlerl fl998t lHughesll200l . Therefore, 
multiple exposures of an X-ray telescope may be necessary 
for some sources. A possibly confusing sour ce of X-ray 
flares are tidal disruptions of mainsequence stars (iRees I 19901 
Cannizzo, Lee, & Goodman 1990; Komossa et al. 2004). De- 
tection of the afterglow of MBH coalescence will help pin- 
point the location and hence redshift of GR sources. The de- 
tailed light curve of the afterglow will probe the structure of 
the accretion disk as it moves toward the black hole, and will 
also shed light on the cosmological assembly of MBHs. 

We thank Aaron Barth, Tsvi Piran, and Tom Prince for valu- 
able discussions. M. M. was supported at Caltech by a post- 
doctoral fellowship from the Sherman Fairchild Foundation. 
E. S. P. was supported in part by NASA ATP grants NAG5- 
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